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Transients  in  the  Vibrational  excitation  of  Cyclobutane 
Using  the  Variable  Encounter  Method'1 

M.  C.  Flowers/  F.  C.  Welters,  B.  D.  Barton1 
and  B.  S.  Rabinovitch 

Department  of  Chemistry  BG- 1 0 
University  of  Washington 
Seattle,  Washington  98195 

Abstract 

The  probability  of  reaction  of  cyclobutane  molecules  in  a 
fixed  time  interval  after  experiencing  a  known  number  of  collisions 
with  a  hot  surface  at  temperatures  between  749k  and  1126K  has  been 
determined  using  the  Variable  Encounter  Method.  Calculations  util¬ 
izing  exponential  or  gaussian  models  for  energy  transfer  enable  the 
average  amounts  of  energy  transferred  for  deactivating  collisions, 
v.\E's,  to  be  estimated.  The  exponential  model  fits  the  experimental 
data  the  best  and,  using  this  model,  <AE'S  is  2430  cm  ^  at  748k  and 
decreases  to  1470  cm  ^  at  1123k. 
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Introduction 

A  new  technique  (the  Variable  Encounter  Method)  has  been  described 
recently  [1]  for  the  study  of  the  relaxation  of  vibrational  energy  of  gaseous 
molecules  on  collision  with  a  hot  surface.  VEM  has  previously  been  applied  to 
the  study  of  cyclopropane  isomerization  [1],  By  this  method,  molecules  are  colli 
sionally  equilibrated  with  a  cold  surface  of  known  temperature  and  then  randomly 
enter  a  hot  reactor  where  they  experience  a  series  of  collisions  with  a  sur¬ 
face  van  encounter).  Molecules  leave  the  reactor  at  random  times,  usual  Is 
without  reaction,  and  are  then  re-equilibrated  with  the  cold  surface  before 
re-entr\  into  the  reactor.  The  average  number  of  collisions  a  molecule 
experiences  with  the  surface  of  the  hot  reactor  on  each  encounter  can  be 
varied  by  changing  the  dimensions  of  the  reactor.  The  fraction  of  reaction 
that  has  occurred  after  a  known  number  of  encounters  can  be  determined.  This 
can  be  used  to  deduce  information  on  the  rates  of  population  of  vibrational 

states  above  the  threshold  energy  for  reaction,  l  ,  i.e.,  in  the  transient 

o 

region . 

In  this  paper  we  describe  the  application  of  VEM  to  the  study  of  the 
vibrational  energy  relaxation  in  cyclobutane.  This  molecule  decomposes  to 
ethylene  and  the  reaction  has  been  well  studied  as  a  homogeneous  thermal 
process  [2-7],  The  Arrhenius  parameters  are  very  similar  to  those  for  cyclo¬ 
propane  so  that  an  interesting  comparison  with  the  earlier  cyclopropane  results 


is  afforded. 
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Experimental 

The  reaction  vessel  consisted  of  a  2  C  fused  quartz  flask  to  which  tubular- 
reactor  fingers  of  the  dimensions  indicated  in  Fiq.  1  were  attached.  The  entrance 
area  of  each  tubular  reactor  was  $  1  of  the  total  area  of  the  reaction  vessel. 

Molecules  randomly  enter  a  reactor  from  the  main  flask.  Stochastic 
calculations,  described  in  ref.  1,  were  used  to  estimate  the  mean  number  of 
collisions,  m  ,  a  molecule  would  suffer  with  the  wall  of  a  particular  tubular- 
reactor  before  it  re-emerged  into  the  flask.  The  values  of  m  were  27.2  for 
reactor  (i),  8.47  for  reactor  (ii)  and  2.6  for  reactor  (iii)  (Fig. 1).  In  the 
case  of  reactor  (i),  molecules  that  collided  with  the  walls  of  the  antechamber 
(A)  but  did  not  then  enter  the  main  part  of  the  reactor  (R)  may  be  assumed  not 
to  contribute  to  reaction.  Molecules  that  left  (B),  collided  with  (A)  and  then 
re-entered  (R)  were  included  in  the  calculation  with  the  assumption  that  the 
temperatures  of  the  walls  of  (A)  and  (P)  were  the  same. 

The  tubular  reactors  were  enclosed  in  stainless  steel  furnace  blocks,  which 
were  heated  by  clam-shell  kanthal -wound  heaters.  The  internal  diameter  of  the 
stainless  steel  block  was  2-3  mm  greater  than  the  outer  diameter  of  the  tubular 
reactors.  Reduced  temperature  gradients  were  obtained  by  filling  the  gap  between 
the  reactor  and  stainless  steel  block  with  fine  ('v  100  mesh)  silica  sand. 

The  top  1  cm  of  reactors  (i)  and  (ii)  were  heated  with  auxiliary  heaters 
made  from  16  gauge  kanthal  wire  set  in  alundum  cement.  Additional  heat  to  the 
top  of  the  stainless  steel  block  for  reactor  (iii)  was  supplied  via  16-gauge 
coi’s  of  kanthal  wire  set  in  cement  in  a  recess  in  the  top  of  the  block.  Maxi¬ 
mum  temperature  gradients  for  reactors  (i)  and  (ii)  were  less  than  ±  25l'C,  with 
75‘T  of  the  area  of  the  reactors  within  t  10°C  or  less.  Temperature  gradients 
for  reactor  (iii)  were  less  than  t  10°C.  Temperatures  were  measured  using 
chromel/alumel  and  platinum/platinum-10‘l  rhodium  thermocouples  cemented  to  the 
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exterior  surface  of  the  reactors.  Temperatures  were  controlled  to  ♦  1°C. 

The  reaction  vessel  was  evacuated  using  a  conventional  glass  diffusion 
pump  to  %  1x10  ^  mm  Hq,  as  measured  by  a  cold  cathode  discharge  gauge.  Cyclo¬ 
butane,  whose  major  impurity  was  0.2  total  of  propene  and  butene-1,  was  thor¬ 
oughly  degassed  before  use.  Prior  to  kinetic  runs  the  reactors  were  "aged"  by 
pyrolysis  of  cyclobutane  at  a  pressure  between  3x10'  and  3x10’  torr  for 
prolonged  periods,  at  a  temperature  equal  to  the  highest  temperature  used  for 
that  reactor.  No  difficulty  was  experienced  with  aging  reactors  (i)  and  (ii). 
Reproducible  rate  constants  for  the  formation  of  ethene,  with  only  trace  quan¬ 
tities  of  propene  and  1-butene  being  observed,  resulted  from  aging  for  48  hours. 
Aging  of  reactor  (iii)  was  more  difficult  and  required  several  days  before  repro¬ 
ducible  rates  were  obtained.  Propene  and  1 -butene  formation  appears  to  be  a 
concurrent  parallel  surface-catalysed  reaction,  in  agreement  with  a  previous 
study  [5].  Reactor  (iii), when  fully  aged,  gave  yields  of  propene  and  1 -butene 
separately  that  amount  to  }  5.  of  ethene  yields. 

The  reaction  was  studied  by  adding  aliquots  of  reactant  to  the  reaction 
vessel  from  a  gas  burette,  pyrolysing  the  reactant  for  measured  times  and  then 
rapidly  freezing  the  reaction  mixture  into  a  trap  containing  a  small  amount  of 
60-80  mesh  Porapak  Q  cooled  to  -196°C.  Ethene  product  was  trapped  quantitatively. 
The  relaxation  time  for  pressure  decrease  during  freeze  down  was  %  10  s.  A  pump- 
down  time  correction  was  added  to  the  pyrolysis  time  used  to  calculate  rate 
constants. 

Analysis  of  the  reaction  products  was  carried  out  by  gas  chromatography 


using  a  6  ft. -1/8"  i.d.  column  of  28‘T  squalane  on  60-80  mesh  firebrick  with 
f.i.d.  detection. 


5 


Results  and  Calculations 

The  rate  of  decomposition  of  cyclobutane  to  ethene  was  investigated  over 

-4 

a  temperature  range  of  749k  to  1126k  at  an  initial  pressure  of  \  1.6x10  ton 
cyclobutane.  At  each  temperature  plots  of  log  (fraction  of  cyclobutane  unreacted) 
ys  time,  which  were  accurately  linear  and  went  through  the  origin,  were  used  to 
determine  apparent  first  order  rate  constants  for  decomposition.  Normally,  the 
extent  of  reaction  at  from  6-8  different  reaction  times  was  determined.  A  typi¬ 
cal  first  order  plot  is  shown  in  Fig.  2.  The  standard  deviation  of  the  experi¬ 
mental  slopes  of  the  first  order  plots  was  less  than  2T .  A  mass  balance  of  better 
than  98  was  obtained. 

A  series  of  runs  were  carried  out  at  1033k  in  reactor  (i),  over  the  range 
-5  -4 

between  9x10  torr  and  5x10  torr  initial  cvclobutane.  No  change  in  the 

measured  first  order  reaction  rate  constant  was  observed.  Runs  in  reactor  (i) 

-4 

in  which  cyclobutane  and  cyclopropane,  at  partial  pressure  of  -v  2x10  torr, 
were  co-pyrolysed  similarly  showed  no  change.  These  mixed  pyrolyses  provided  a 
very  useful  comparison  with  previous  work  on  cyclopropane  [1]  and  gave  comparable 
results  [8], 

From  a  knowledge  of  the  dimensions  of  the  entrance  orifice  to  each  reactor, 
the  value  of  m  ,  the  molecular  weight  of  cyclobutane,  and  the  mean  temperature 
of  the  main  reactor  vessel,  the  apparent  first-order  rate  constants  were  converted 
into  mean  probabilities  for  reaction  per  collision  with  the  reactor  wall,  ^(m). 
Values  of  Pc (m)  are  given  in  Table  I  and  values  of  log  P  (m)  vs  temperature  are 
plotted  in  Fig.  3. 

In  order  to  deduce  information  on  energy  transfer  between  cyclobutane 
molecules  and  the  hot  surface  from  the  experimental  data,  several  models  for 
the  energy  transfer  probability  distribution  function  were  explored.  A  com¬ 
puter  simulation  of  the  collisional  processes  was  carried  out  and  comparisons 
made  between  calculation  and  experiment. 
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The  computer  modelling  procedure  has  been  described  in  detail  in  ref.  1. 

Two  different  ("flat")  models  for  the  probability  of  a  given  down  transition, 

PA£  were  used:  Model  A  (exponential);  P  ^  =  A  exp  (-AE/<AE>)  for  0  <  AE  c  9000 
cm"1,  P„£  =  0  for  AE  >  9000  cm"1;  Model  B  (Gaussian);  P^  =  A'  exp  {-(AE-  AEmp)^/2c^} 
for  0  $  AE  $  9000  cm  1 ,  P.^  =  0  for  AE  >  9000  cm  1 .  Here,  A  and  A'  are  normali¬ 
sation  constants;  <AE>,  AEmp  (mp  signifies  most  probable)  and  a  are  parameters 
of  the  model;  AE ,  AEmp  were  constant  ("flat")  models,  and  a  was  set  equal  to 
0.7  AEmp.  The  truncation  AE  <  9000  cm  1  is  a  practical  computational  feature  to 
limit  the  transition  probability  matrix  to  tractable  dimensions.  Detailed  balance 
and  completeness  were  maintained.  It  should  be  noted  that  because  of  truncation, 
the  effective  average  down  transition  size,  called  <AE’>,  is  not  equal  to 
<AE>  ,  except  when  the  latter  quantity  is  small;  <AE ' >  may  be  significantly 

different  from  <AE>  and  AE  . 

mp 

The  microscopic  rate  constants,  k^ ,  for  decomposition  of  cyclobutane  mole¬ 
cules  with  a  specific  energy  content,  which  are  required  in  the  computer  simulation, 
are  calculated  from  RRKM  theory.  The  reaction  coordinate  was  assumed  to  be  a 
C-C  stretch  of  the  cyclobutane  molecule  and  the  frequencies  used  for  the  transition 

state  (App.  A)  were  chosen  to  give  agreement  with  the  high  pressure  rate  constant 

15  6 

at  772K  as  calculated  from  the  Arrhenius  equation  [3],  k^  =  10  ‘  exp 

(-62500/1. 987T)  sec"1. 

The  results  of  the  computer  simulation  are  included  in  Fig.  3  for  both 
models.  In  each  case  a  value  of  <AE ' >  was  found  that  produced  a  fit  to  the 
experimental  data  for  reactor  (iii)  (m  =  2.6).  The  curves  for  reactors  (i) 
and  ( i i )  were  then  calculated  using  the  same  values  for  <AE’>.  This  approach 
was  adopted  since  the  value  of  Pc(m)  for  the  reactor  with  the  smallest  mean 
collision  number  is  the  most  sensitive  to  changes  in  <AE>. 
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Calculations  using  "linear"  gaussian  and  exponential  models  in  which  IE 

mp 

or  <AE^  (and  thus  <AE 1  >)were  allowed  to  vary  with  initial  energy  according  to 
the  relationship  <AES  (or  AE  )  =  a  +  6E ,  where  a  and  6  were  adjustable  parameters, 
were  also  carried  out.  These  calculations  gave  results  in  poorer  agreement  with 
the  experimental  data  for  reactors  (i)  and  (ii),  particularly  at  the  highest 
temperatures.  The  calculated  values  of  Pc(m)  for  higher  m  are  somewhat 
larger  than  those  calculated  for  the  flat  models  in  which  <AE'>  does  not  vary 
with  initial  energy. 


Di scusss ion 


Conventional  studies  of  the  thermal  decomposition  of  cyclobutane  [2-7] 
indicate  that  the  reaction  is  a  homogeneous,  unimolecular  reaction,  although 
at  the  lowest  temperatures  investigated  earlier,  the  reaction  rate  did  not 
decline  as  rapidly  as  would  have  been  expected  for  a  unimolecular  reaction  in 
the  fall-off  region  [5,6].  This  behavior  was  probably  due  to  the  increasing 
importance  of  surface  activation  at  low  pressure.  Indeed,  surface  activation 
has  been  shown  to  be  important  in  the  thermal  decomposition  of  methyl  cyclo¬ 
butane,  but  with  no  evidence  for  surface  catalysis  [9]. 

In  the  present  study,  we  believe  also  that  the  only  process  of  significance 
occurring  at  the  wall  is  energy  transfer.  This  allows  an  estimate  to  be  made 
of  the  amounts  of  energy  transferred  on  average  per  collision. 

The  flat  exponential  model  (model  A)  produces  the  best  agreement  with  our 
experimental  data  over  the  entire  temperature  range  studied.  This  contrasts 
with  a  previous  conclusion  [1]  from  VEM  studies  with  cyclopropane  that  flat 
gaussian,  or  a  linear  exponential,  gave  the  best  fit;  however,  the  change  in 
model,  if  real,  is  in  line  with  the  somewhat  decreased  stepsize  measured  in  the 
present  study  (e.g.,  <AE '  >  -  2000  cm’^  vs  <AE'>  =  2990  cm  "*  for  cyclopropane,  at 
950K,  both  with  a  gaussian  model).  The  uncertainties  are  such  that  it  is  prob¬ 
ably  not  possible  to  decide  definitively  what  the  exact  form  of  the  correct  model 
is. and  at  the  lowest  temperatures  the  distinction  between  the  models  is  small. 

Average  amounts  of  energy  transferred  per  deactivating  collision  (<AE’>) 
do  not  change  dramatically  between  the  two  models  as  can  be  seen  from  Fig.  3. 

The  actual  values  of  <AE’>  are  dependent  on  the  absolute  values  of  P  (m)  and  are 
subject  to  some  experimental  uncertainty. 

Plotted  in  Fig.  4,  as  a  function  of  energy,  are  the  ratios  of  the  transient 
population  numbers  of  molecules,  N(E),  in  energy  state  E,  to  the  numbers, 
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N(E)eq  ,  that  would  be  present  if  the  molecules  were  in  thermal  equilibrium 
with  the  reactor  surface  with  no  reaction  occurring.  The  ratios  are  shown  after 
various  numbers  of  collisions  with  the  surface.  It  can  be  seen  that  at  the  low¬ 
est  temperature,  progression  to  the  steady  state  occurs  much  more  rapidly  than 
at  the  highest  temperature.  The  marked  depletion  of  molecules  at  energies  above 
the  critical  energy  for  reaction  is  the  consequence  of  the  chemical  reaction. 

The  curves  in  Fig.  4  were  calculated  for  the  exponential  model;  however,  the 
qualitative  picture  is  the  same  using  the  gaussian  model. 

Incubation  times,  in  units  of  number  of  collisions,  are  defined  as  the 
difference,  t.^  =  tfp  -  t*p  ,  between  the  mean  first  passage  time  and  the  steady 
state  mean  first  passage  time;  xinc  is  9.4  and  7.4  at  748K,  and  18.3  and  16.1 
at  1123K,  calculated  on  the  exponential  and  gaussian  models,  respectively.  For 
cyclopropane  [1],  incubation  times  were  found  to  be  a  little  smaller. 

Plots  of  the  probability  of  reaction  after  the  nth  wall  collision  for  models 
A  and  B  are  shown  in  Fig.  5.  At  low  values  of  n  ,  the  exponential  model  gives 
higher  reaction  probabilities  per  collision  than  the  gaussian  model,  but  this  is 
reversed  at  higher  values  of  n  .  Reactors  with  small  values  of  m  ,  of  course, 
weight  the  small  collision  numbers  more  highly  than  large,  and  it  can  thus  be 
seen  why  the  exponential  model  requires  the  smaller  value  of  <AE'>  to  fit  the 
data  for  m  =  2.6;  but  this  method  of  fit  then  produces  a  higher  mean  probability 
of  reaction  per  collision  in  reactors  with  m  =  8.47  and  27.2,  on  the  basis  of 
the  gaussian  model,  than  is  found  for  the  exponential  model. 

There  appears  to  have  been  no  previous  studies  for  which  weak  collider 
calculations  have  been  made  for  cyclobutane.  Values  of  <AE>  have  not  been 
reported.  Estimated  values  for  the  collisional  efficiency  B  of  various  gases 
for  energy  transfer  for  gas-gas  collisions  relative  to  cyclobutane,  have  been 
reported  around  722K,  as  based  on  the  shifts  in  k/k^  curves  on  addition  of  inert 
gases  [4].  The  values  of  6  are  approximately  the  same  as  those  found  for  cyclo¬ 
propane  [10]  using  the  same  procedure. 
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In  thermal  studies  it  is  well  established  that  complex  (large)  bath 
molecules  and,  of  course,  the  parent  substrate  tend  to  behave  as  strong 
colliders  at  low  temperature  [11].  Recent  chemical  activation  work  at  300*., 
with  cyclobutane  as  the  bath  gas,  has  tended  to  favor  high  values  of  <AE> 

(=  10  kcal  mole'1 ) for  the  collisional  deactivation  of  highly  vibrationally 
excited  methylcyclobutane  [12]. 

The  important  conclusion  reached  by  this,  and  by  the  previous  study  on 
cyclopropane  [1],  is  that  the  surface  behaves  as  an  increasingly  strong  colli¬ 
der  as  the  surface  temperature  is  decreased.  In  general,  at  temperatures  in 
the  vicinity  of  those  used  in  conventional  kinetic  studies  of  collisional 
efficiency,  i.e.,  <  BOOK,  wall  collisions  are  more  effective  for  energy  transfer 
than  are  gas-gas  collisions. 
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Table  1.  Probabilities  of  reaction  of  cyclobutane  per  wall  collision  (P^(m) 
in  reactors  (i),  (ii),  and  (iii) 


Reactor  (i) 

Temp(K) 

1080 

1034 

960 

890 

843 

790 

749 

(m  =  27.2) 

Pc(m)xl07 

745 

331 

75.4 

13.8 

3.53 

0.506 

0.0783 

Reactor  (ii) 

Temp(K) 

1075 

1033 

958 

894 

8^0 

792 

(m  =  8.47) 

Pc(m)xl07 

283 

148 

31.4 

6.57 

1.37 

0.206 

Reactor  (iii) 

Temp(K) 

1126 

1072 

994 

907 

823 

(m  =  2.6) 

P  (in) xl 07 

19.2 

9.54 

2.55 

0.413 

0.0514 
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Appendix:  Parameters  used  in  RRKM  calculation  for  the  decomposition  of 

cyclobutane  to  ethylene. 

Cyclobutane  frequencies  (cm’^ )  [13]:  2987,  2975,  2952(2),  2895,  2893, 

2887(2),  1447(2),  1443(2),  1260,  1257(3),  1223(2),  1222,  1219,  1001(2), 
926,  898(2),  7492(2),  741,  627,  197. 

Transition  state  frequencies  (cm  ^):2980(8),  1447(2),  1443(2),  1260(2), 
1250(2),  1100(2),  1000,  898,  850,  740(2),  600(2),  380,  350(2),  150 
Reaction  path  degeneracy  =  4. 


_ * 


_ ■ 
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Figure  Captions 

Fig.  1  Detail  of  tubular  reactor  fingers  attached  to  2  «  fused  quartz 
flask,  (all  dimensions  are  in  mm) 

Fig.  2  First  order  plot  for  the  decomposition  of  cyclobutane  to  ethylene 
at  994K  for  the  m  =  2.6  reactor. 

Fig.  3  Plots  of  experimental  values  of  log^  Pc(m)  vs  temperature,  and  of 
calculated  curves  for  models  A  and  B  for  each  reactor.  Experimental 
data  (each  point  is  the  result  of  at  least  6  separate  runs): 

O  .  m  =  27.2;  ,  m  =  8. 47 ;  Q  ,  m  =  2.6.  Solid  curves  are  the 

calculated  curves  using  model  A  (exponential)  and  the  dashed  curves 
are  those  for  model  B  (gaussian).  Values  of  <AE‘>  used  are  indicated 
at  various  temepratures.  The  parenthetic  E  and  G  designate  the  expon¬ 
ential  and  gaussian  models,  respectively. 

Fig.  4  Ratio  of  population  of  molecules  relative  to  the  equilibrium  population 
(at  the  temperature  of  the  reactor  wall)  as  a  function  of  energy, 
using  model  A.  The  number  of  previous  collisions  experienced  by  the 
molecules  is  shown  alongside  the  curves;  (a)  reactor  wall  =  748K, 
<AE’S  =  2430  cm  \  (b)  reactor  wall  =  1123K,  <AE’>  =  1470  cm'1. 

Fig.  5  Histogram  of  probabilities  of  reaction  after  collision  n  ,  and  before 
the  next  collision,  as  a  function  of  n  for  model  A  (solid)  curve  and 
model  B  (dashed)  curve;  (a)  748K;  model  A  <AE’>  =  2430  cm’1;  model  B 
<AE ' >  «  2830  cm'1  (b)  1123K;  model  A  <AE'>  «  1470  cm"1, 
model  B  <AE'>  -  1860  cm'1 
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